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We study models of fermion mass matrices based on a flavor- and generation-dependent 
string-motivated U(l)^ gauge symmetry and report two new classes of solutions to the 
requisite consistency conditions. In particular, we propose that the fundamental reason 
underlying the striking feature mj, m T «m t is that all of the elements of the down-quark 
and charged lepton effective Yukawa matrices actually arise from higher- dimension operators, 
suppressed by inverse powers of the Plank mass. An explicit model embodying this idea is 
H ' constructed. 
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The pattern of fermion masses and mixing remains one of the most important mysteries 
in particle physics. The successful standard model (SM) can accomodate but not explain 
this pattern. A satisfactory understanding would require that one have an experimentally 
confirmed theory explaining the related electroweak symmetry breaking (EWSB), and one 
does not have this at present. Nevertheless, one may proceed by exploring plausible models. 
The fermion mass spectrum has several striking features: (i) within each charge sector, the 
masses increase with generation by large factors: m u « m c « m t , m d « m s « mi,, 
and m e « m M << m T ; (ii) if one assumes that all of these masses arise from conventional, 
dimension-4 Yukawa couplings, the associated Yukawa couplings for all of these fermions 
except the top quark are all much smaller than a typical small coupling like e = yAna ~ 0.3, 
without explanation; (iii) related to (ii), even if one restricts to the third generation, the 
masses are still quite different: m T and m& are both « m t . A related feature is that (iv) 
the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix is near to the identity. If 
one assumes certain simple forms for Yukawa matrices with zeros in various entries, it is 
possible to explain property (iv) as a consequence of (i), because the quark mixing angles 
are functions of (square roots of) small mass ratios like md/m s , m u /m c , etc. [p], ||. However, 
in such an approach, the fermion masses are used as inputs, and properties (i), (ii), and 
(iii) are not explained.^] Indeed, although there has been recent progress in explaining (i) 
and (iv) via contributions of higher-dimension operators at a high mass scale near to the 
scale of quantum gravity [[|, such efforts have not addressed what, to us, seems an equally 
remarkable feature, viz., (iii). The attractive idea of radiative electroweak symmetry breaking 
in a supersymmetric generalization of the standard model j|] depends on the existence of at 
least one quark which has a mass comparable to the EWSB scale, but it cannot explain why 
this was the top quark instead of the bottom quark (or, indeed why both m t and m;, are not 
~ the EWSB scale), and hence it cannot explain (iii) or the full extent of (ii). 

In this paper, we shall explore an appealing class of models of fermion mass matrices 
which has the potential to explain all of the properties (i)-(iv). We shall present a particular 
model which, we believe, is the first to offer a possible fundamental explanation of property 
(iii). The explanation is that in the down-quark and charged lepton sectors, the masses 
of not just the first two generations, but of all generations arise from higher-dimension 
operators suppressed by powers of a small mass ratio, e oc v/Mp, where v is the breaking 
scale of a flavor- and generation-dependent, string-motivated U(1) J 4 symmetry, and Mp = 
(87rG JV )- 1 /2 = 2.44 x 10 18 GeV is the (reduced) Planck mass. 

To set our work in context, we note that since the original success of the standard model, 
there has been a growing appreciation that its renormalizability and, in particular, the 
absence of higher- dimension operators, may well be the consequence of a large logarithmic 
interval in energy between the electroweak scale and a higher scale where there is new 
physics (see, e.g., Ref. 0). Indeed, there are specific reasons for expecting such operators 
at this high scale: the only known way to stabilize the hierarchy v EW << Mp (where 
v EW = 246 GeV is the EWSB scale) is via a supersymmetric generalization of the SM. In 
turn, global supersymmetry is naturally embedded in supergravity, which one also finds as 

1 Indeed, the model of Ref. [9 is excluded by property (iii), since in that model, \V c b\ — |\/ m s/ m b — 
'y 'm c j l mt\ (masses normalized at same scale (j,); because m t » m c , the second term is too small to 
significantly cancel the first term -J m s /mb ~ 0.17 and fit the experimental value, \V c b\ = 0.038 ± 0.005. 
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the low-energy limit of the main candidate for quantum gravity, string theory. But d = 4 
supergravity is nonrenormalizable. Indeed, explicit calculations of the pointlike limit of string 
theories for energies E « M str (where M str = 2(c/)~ 1 / 2 = gMp) yield supergravity as a low- 
energy effective field theory with infinite towers of higher- dimension operators with coefficient 
functions proportional to inverse powers of M str and powers of the compactification scale (the 
latter may be non-explicit in four-dimension string formulations). One must therefore take 
account of higher-dimension operators when analyzing terms which contribute to fermion 
mass matrices. This was, indeed, already realized long ago 0, although detailed studies 
have only been carried out recently. At first, one might consider these higher-dimension 
operators to be an unfortunate, if inevitable, complication in the theory. However, they 
may well play a very important role in the area of fermion masses. Specifically, via vacuum 
expectation values (vev's) of the scalar components of certain chiral superfields, which we 
shall denote generically as v, these higher-dimension operators can yield contributions to 
effective dimension-4 Yukawa interactions which are suppressed by powers of the ratio e ~ 
v/M str . This idea has already been used for a possible explanation of properties (i) and (iv) 
||; here we shall extend this work with new solutions of the consistency conditions and take 
a step further, to use the small ratio e to explain (ii) and (iii). A preliminary report of some 
of our results was given in Ref . |7J . Related results in a somewhat different direction (having 
nib and m T arise from dimension-4 operators, as in Ref. ||) were presented in Ref. ||. 

We shall work within the context of a supergravity theory which reduces at low energies 
to the minimal supersymmetric standard model (MSSM). We consider a theory where there 
is a flavor and generational symmetry group Gp which restricts the forms of the terms in 
the action. In particular, this symmetry forbids certain cubic superfield couplings which give 
rise to Yukawa interactions. Since this happens at a scale not too far from that of quantum 
gravity, and since, in general, global symmetries are broken by quantum gravity || even at 
the semi-classical level, one is motivated to make Gp a gauge symmetry. One then faces two 
questions: (a) is there a natural origin for Gf in the presumed underlying string theory? 
and (b) is there a natural way to explain why the scale of the breaking of Gf is such that the 
ratio v I M str has the value that it must to fit the observed forms of fermion mass matrices? 
A possible affirmative answer to both of these questions is provided by a gauged symmetry 
Gp = U(1)a which is, at the field theory level, apparently anomalous, but whose anomaly is 
cancelled by a Green-Schwarz mechanism |f[0|| . Such U(1)a gauge symmetries are known to 
arise is various string models and, moreover, they are broken at a calculable scale v given by 
v 2 ~ M 2 tr /(192n 2 ) |TJ, so that e ~ (Sttv^) -1 = 0.023, a value which is in the right general 
range to explain fermion mass hierarchies 0. Of course, such a 11(1)^ symmetry does not 
mix up-type and down-type quark superfields, or mix these with lepton superfields and is 
thus quite different from flavor and generational symmetries which comprise extensions of 
grand unified groups. 

We shall denote the (left-handed) SM matter chiral superfields as Qi, u% d\, Li, and ef, 
where i = 1,2,3 is the generation, with {ui} = {u, c, t}, {di} = {d, s,b}, etc. Under the 
flavor- and generation-dependent U(l)x, these carry the charges q Q , q u c, q d c, q L , and q e c. 
The Y — 1, — 1 Higgs chiral superfields are denoted Hi and H 2 and have U(l)^ charges q 
and q H . We shall also assume that the theory is invariant under the usual R parity. The 
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cubic superfield terms in the globally supersymmetric superpotential are then given by 

Wcwc = Q i Y^ ) u c J H 2 + QiYWq H x + L i Y®e c j H 1 (1) 

We shall assume that the gauge symmetry at energies > v is G = G$m x U(1)a, where G$m 
is the SM gauge group SU(3) x SU(2) x U(l)y. The theory will also contain certain chiral 
superfields which are SM-singlets but transform under the U(l)^. A discussion of constraints 
on these to avoid destabilization of the gauge hierarchy is given in Ref. [T2 . 



To proceed, we recall how the cancellation of an apparent field-theoretic anomaly in a 
gauged U(l) symmetry works when this U(l) is, in fact, non-anomalous in the full string 
theory ( Green- Schwarz mechanism) |T(J and the related (Dine-Seiberg-Witten, DSW) mech- 
anism whereby the U(l) is broken [|ll||. Here, we shall denote chiral superfields generically 
by $. Standard d = 4 supergravity is described by two functions. The first is the generalized 
Kahler potential 

= K($,$) +1bW($) + \nW(&) (2) 

is a hermitian function of the Kahler potential K and the superpotential W, the latter being a 
holomorphic function of the chiral superfields. The second is the gauge kinetic normalization 
(matrix) function / a fc($), where a,b = 1, ..,Ng, the number of gauge bosons. As indicated, 
f a b is a holomorphic function of the chiral superfields. We shall use G, K, and W to denote 
both superfield quantities and also their scalar components; however we shall distinguish 
between chiral superfields and their scalar components by type case (e.g. <ft = m = 
The bosonic part of the standard supergravity lagrangian is pT 



-^C B = GijD^g^D^ + ^R- hief ab F« u F b ^ - hmf ab F« u F b ^ - V, (3) 

where the auxiliary field contributions to the potential are given by V = V + D, with 

V = e G {G i G ] G fj -3) (4) 

and 

V = Re(f^)D a D b , (5) 

where D a is the D-type auxiliary field associated with generators T a of the gauge group, 
normalized according to Tr[T a ,Tb] = (l/2)5 a &. Discussions of one-loop corrections are in 



Ref. |15| and references therein. The scalar fields and their conjugates are denoted as above 
by (f) % and <p % = (jf. The Kahler metric Gq = Kq = d^d^ has inverse G l K Finally, the 
derivatives D^, which are gauge and general coordinate invariant, are normalized so that 

D tt = d lt + z a iA%r a . 

The function f ab determines the gauge coupling constants of the theory and also plays 
an important role in the cancellation of the field theoretic anomaly due to the FF term in 
the Lagrangian. In string models, at tree level, f a b is given by 

fab = k a s5 abj (6) 

where the levels k a of the Kac-Moody algebras on the worldsheet depend on the gauge 



factor group G a [16|. We denote ki, i = 1, 2, 3 and kj± respectively as the Kac-Moody levels 
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corresponding to the factor groups U(l)y, SU(2), and SU(3) of Gsm and the flavor symmetry 
U(1) A . The dilaton vev determines the gauge couplings, according to l/g% = k a < Re(s) >. 
At one- and higher-loop levels, f a b acquires a dependence also on the moduli fields ]17|, [18], [TJ§ 
from chirally anomalous triangle diagrams in the effective field theory as well as string 
threshold corrections generically needed to cancel these anomalous terms, as is required 
by the underlying string theory. However, there is a large class of models in which the 
field theoretic chiral anomalies can be cancelled by a "universal" Green-Schwarz mechanism 
involving the dilaton, similar to the gauge U(l) anomaly cancellation, in which the one-loop 
corrected /„& is such that Re(f a b) (which is what controls the unification conditions for the 
gauge couplings) has the form Re(fab) = k a 8 a b(Re(s) + f(T, T)), where T denotes the moduli 
|T7| , |1~8| , [I9| . We shall restrict our attention to this class of models. (There are also small 
logarithmic corrections from the running of the gauge couplings between M stT and Mx, the 
scale of gauge coupling unification; as in [Q, we shall neglect these here since M x is not 
« M str .) 

If the effective theory contains no gauge anomalies, then (|3]) is gauge-invariant, with the 
dilaton being a gauge singlet. In the case of interest here, however, where G contains a U(1)a 
factor which has field theoretic anomalies, then one-loop corrections from light (< Mpi) 
fermions give an anomalous correction C anom to ([3|). Under a U(l)^ transformation, C anom 
transforms by 

5C anom ex caFFF* + £ CaFrP^ + (7) 

a 

where 

c A = (l/3)Tr(T A T A T A ) (8) 

c a = Tr (T A T a T a ), (9) 

and a runs over the gauge group factors which have no field theoretic anomalies (in the 
present case, the factors in Gsm)- The ellipses in (5) indicate additional terms which we 
shall discuss below. 

The theory is invariant under constant rescalings — > aA%, which amounts to a simul- 
taneous redefinition of k a and T a : k a — > a 2 k a and T a — > aT a . (Our normalization of the 
generators of Gsm is the standard one.) For the U{1) A group the above rescaling could be 
used to set k A to 1 at the expense of modifying the U{1) A charges. We shall instead use this 
rescaling to fix the U(1) A charge of one of the SM-singlet particles. 

Due to the coupling of Im(s) to FF in (||), this variation of C anom can be cancelled by 
assigning a U{1) A gauge variation to Im(s) |10|, |TT|, but only if 

c A : c a :: k A : k a for all a. (10) 

Since s is no longer a gauge singlet, its kinetic terms must be made gauge invariant. At 
the superfield level, this means that the Kahler function K must include a coupling between 
the chiral superfield S and V A , the real superfield that contains the U{1) A gauge multiplet. 



Specifically, the tree level Kahler function for S is modified as WY 
-\n(S + S) -> -\n(S + S + cV A ) 
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where c is a constant determined by gauge invariance and is related to the coefficients ®. 
The second term in the last line of ( |iT| ) produces a term in Lb which is linear in the auxiliary 
D field of Va- Once Re(s) acquires a nonzero vev, this causes some of the fields which are 
charged under the U(1)a to acquire vev's, thereby spontaneously breaking 11(1)^- 

In a given string model, ( |T0D is automatically satisfied, due to the absence of gauge 
anomalies in the underlying theory. The absence of anomalous terms that mix the different 
gauge group factors requires that 

Tr(T A T A T a ) = for all a. (12) 

This ensures that 5C anom does not contain any FF terms that cannot be cancelled by a gauge 
variation of s. Furthermore, mixed gauge and gravitational anomalies must be absent. In 
addition to terms explicitly displayed in ([□]), the U(1)a gauge variation of C anom contains 
a term proportional to Ti(TA)R flu Rf 1 u- This term can be cancelled because the effective La- 
grangian contains a nonstandard (higher derivative) coupling /cgIhi(s).R AU/ 'R pv (which defines 
ka)- The cancellation requires a relation between the coefficients c a ,CA and Tr(T^) because 
the same U(1)a transformation of s must cancel all anomalies. This last relation is difficult 
to use phenomenologically because (i) Tr(T^) depends on SM singlets about which we have 
very little experimental information, and (ii) because it depends on fc^ and ka which are not 
measured at low energies. 

We proceed to investigate such models to obtain phenomenologically acceptable fermion 
mass matrices. Within the above theoretical framework, we shall make the following specific 
assumptions: 



The low energy theory near the electroweak scale is the MSSM with phenomenologically 
viable soft supersymmetry breaking and supersymmetric mass terms. This assumption 
precludes SM nonsinglets which carry U(1)a charge from getting large vev's due to the 
U(1) A D-tenn. 

Supersymmetry is spontaneously broken in a hidden sector with 
m 3/2 = < e G > ~ rriyy and < V > ~ 0; i.e., 0{M P ) and 0(Mpmy 2 ) contributions 
to < V > cancel. This ensures that the soft breaking terms are naturally of 0(mw)- 
In addition, this assumption means that the dominant contributions to the effective 
Yukawa couplings are only from superpotential terms, as in the MSSM, and not from 
K, which can also potentially contribute to the effective Yukawa couplings in a super- 
gravity model |f20 . 



3. The SM gauge couplings unify in the canonical way, i.e., g 3 2 : g 2 2 : g Y 2 = 1:1: 5/3, 
or equivalently, k$ : hi : h\ — 1 : 1 : 5/3. 

4. SM singlet fields \ which get a vev with < Da > = carry U(1)a charges with the same 
sign. In the class of models that we consider, < x > ~ 10 16 GeV, which is >> %. 
This ensures that < W > (and hence < e G >) can be kept small without resorting to 
accidental near-cancellation between two large vev's in < W > or introducing extra 
symmetry. For example, if the U(1)a charge of % is 1 and the U(1) A charge of another 
SM singlet x' is — b, b = 2,3, then gauge invariance allows W to contain terms of 



5 



/ 


<5 




d c 


L 


e c 




9 Q +«i 


9u c + A 


9d c + 7i 


q L + 01 


g e = + &i 


9/ 2 


9q +"2 




9> + 72 


9 L +«2 


g e c + 6 2 


9/3 


9 Q - ("l + "2) 


9«<= - (A + AO 


9> - (7i + 72) 


9z, - (oi + aa) 





Table 1: U(1)a charges of SM matter chiral superfields f{. 



the form x b X* /Mp . If both x an d x' S e ^ vev's ~ Mp/100 then the contribution of 
this term to < W > is phenomenologically too large, unless b is sufficiently large. 

5. To reduce the number of parameters, it is convenient to assume that the effective 
Yukawa matrices are symmetric at Mj. 

The hierarchical structure of the low-energy fermion masses can be produced by Yukawa 
couplings which are either hierarchical or democratic. In the latter case, all entries of a 
Yukawa matrix are the same to leading order, while subleading corrections to these differ so 
that the eigenvalues of the matrix have the desired hierarchical structure. A U(1)a symmetry 
can potentially explain hierarchical Yukawa matrices but not (by itself) democratic matrices. 
In contrast, a hierarchical structure can be produced by assigning U(1)a charges to the 
operators in (|l|) in a manner such that for different i,j, these couple to different requisite 
powers of certain SM singlet chiral superfields. 

We next determine which U(1) J 4 charge assignments for the fields satisfy the various 
constraints and lead to experimentally viable fermion Yukawa matrices. There are N q = 
5Nq + 2 = 17 U(1)a charges in the MSSM, where Nq = 3 is the number of matter 
generations. In addition, there are iV s U(1) J 4 charges for SM-singlet chiral superfields. 
The conditions (0) and (|P2"D , together with the assumption of symmetric mass matrices, 
reduce these 17 parameters to 8. To show this, we first list the U(1)a charges of the SM 
matter fields as in Table [1], where q~f = (1/3) J2i=i 9/j is the generational average charge for 
the / chiral superfield. For the Higgs chiral superfields, we put q{Hj) = q H , i = 1,2. To 
ensure symmetric Yukawa matrices, we require that the U(1)a charges of the chiral superfield 
bilinears satisfy 

q(QiU°) = qiQri) (13) 
q(Qid°) = q(Q j( %) (14) 
q(Ue<) = q(Lrf) (15) 
In all, there are six independent contraints in these equations. The solutions are 

«i = A = 7i > &i = h , for i = 1, 2 (16) 

Because the anomalies q in eq. ([|) are linear in the U(1)a charges, it follows that, regarding 
the contributions from the matter fields, the q only depend on the generational average 
U(1)a charges. These anomalies are 

U(1) 2 Y U(1) A : (l/4)[Ar G {(2/3)g Q + (16/3)g uC + (4/3)g d c + 2g i +4g eC } + 2g Hi +2g H2 ] (17) 
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SU(2) 2 U(1) A : c 2 = (1/2) [iV G (3g Q + q J + g Hi + g Ha j (18) 
Sf/(3) 2 t/(l) A : c 3 = (N G /2)(2q Q + q uC + q dc ) (19) 

where we show the general N G dependence but take N G = 3 here. The anomaly conditions 
(|10|) yield two linearly independent constraints on the 7 parameters q Q , q u c, q d c, q L , q e c, q 
and q H2 , e.g., c 2 = c 3 and Ci = (5/3)c 2 . These can be solved in terms of the 5 quantities x, 
y, z, v, and w according to 

q Q = v+x, q u c = 2v+x, q d c = w+y, q L = y, q e c = x, q H± = 3(v+w)+z, q H2 = -z 

(20) 

The U(l)yU(l) A anomaly (|2T|) is quadratic in the U(1)a charges, and hence, in general, 
cannot be written just as a function of the generational averages of the matter field charges. 
However, given (P~5[)- , this anomaly also depends only on these averages; requiring that 
it vanish gives 

U(1) Y U(1) 2 A : = c YAA = N G (f Q - 2q 2 uC + q% -f L + q%) + < - < (21) 



In terms of the 5 quantities in eq. (0), eq. (|21| ) is 

= 2{w 2 - v 2 ) + 3v(w - x) + vz + w(y + z) (22) 



We find the following three families of solutions to (22), which are thus solutions to the total 
set of anomaly constraints (these solutions are independent of N s ): 



qq in? Qd° ql q e ° q H2 q Hl 

xx y y x z —z 

x x \{y — z) y x — \{?>y + z) —z 

x + v x + 2v y + w y x 3(v + w) + z —z 



(23) 



The first two correspond to v = and describe two distinct 3-parameter families of solutions. 
The last exists for v ^ and describes a 4-parameter family of solutions with x solved for in 
( p2]) . The first solution in (^) was already given in Ref. 0]; the other two were not mentioned 
there and are a new result in the present work. The four parameters a\, a 2 , cii, o 2 , together 
with the unknown parameters in (|23| ) yield all allowed 11(1)^ SM charges consistent with 
our assumptions. We find that the constraints are very restrictive, as will be seen. 

In order to account for the important feature (ii) that m t is comparable to the EWSB 
scale v EW , one chooses the source of m t to be a renormalizable, dimension-4 Yukawa coupling, 
as in the SM. This requires the 11(1)^ charge of Q? J u^H 2 to be zero, i.e., 

q Q + q~u° + q H2 = 2(«i + « 2 ), (24) 

Let us denote the 11(1)^ charge of (Q 3 d^Hi) as 9. Then 

q Q + qd° + q Hl = 2(ai + a 2 ) + 9. (25) 
Under this assumption, the U(1)a charges of (QiU C jH 2 ) are given by the matrix 

4a i + 2a 2 3«i + 3a 2 

3«i + 3a 2 4a 2 + 2a\ 2a 2 + ot\ \ (26) 
2«x + a 2 2a 2 + ot\ 
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while 

qiQid'Hx) = qiQiU^) + 9. (27) 

If one took 9 = 0, then the b quark mass would arise from a renormalizable cubic superfield 
operator, and one would not have any fundamental explanation of property (iii). The origin 
of the large mass ratio m t /mb, rather than being explained naturally, would have to be 
pushed into a similarly large value of tan/? = v 2 jv\, viz., tan/? ~ m t /m{,. Instead, we take 
9^0, implementing our explanation of properties (ii) and (iii), since then arises from 
higher-dimension operators, and m& « m t follows naturally. 

We proceed to construct an explicit model for fermion mass matrices embodying these 
ideas. We assume that there are N s = 2 chiral superfields which are SM singlets (and 
U(1)a nonsinglets), X an d x' with unequal U{1)a charges q x and q x '. For convenience, let 
us normalize q x = 1 (this can be done by rescaling the 11(1)^ coupling) and define qy = 9'. 
The DSW breaking of the 11(1)^ yields values for < x > / M P ~ < x' > / M P (denoted e 
above) which are ~ 0(A 2 ), where A = \V US \ ~ 0.22 is a measure of the hierarchical structure 
of the CKM matrix. First, consider the up- and down-quark masses. We have 

W 3 y?j{x, x')H 2 Q l u c 3 + yffix, xlH.Qd'j, (28) 

so that =< yfj > and Y£ =< yf^ >. Motivated by string theory considerations, we 
shall consider only functions y which do not contain fractional or negative powers of x an d 
x' ■ In order to obtain a viable form of the quark mass matrices, we require Y 22 ~ A 4 , i.e. 
(3^2 ~ X 2 , X' 2 or x'x) or (K> M 3 ~ A 2 , i.e. ^23 ~ X or x')- 
We consider three cases: 

(1) 3^22 ~ XX' ■ U(1)a charge conservation (CC) then implies 3^23 = ^32 7^ X>x'- 

(2) y& ~ x- Then CC y u 22 ~ X 2 ■ 

(3) y& ~ X '- Then CC ^2 ~ x' 2 - 

We have studied all of these cases and find in case (2) an assignment of U(l)^ charges 
which yields effective Yukawa matrices (i.e., the matrices which enter in effective dimension-4 
Yukawa terms, all of which, except for Y33, actually arise from higher-dimension operators) 
which give an acceptable pattern of fermion masses at the electroweak scale. At M x , these 
are close to the simple forms 






A 6 





A 6 





A 2 





A 2 


1 



oc A 6 A 2 , Y (d) 



/ A 4 

A 4 A 3 A 3 I , (29) 
V A 3 1 



where the actual entries in the positions given by zeros in ( |2"9"| ) need not be, and are not in 
general, exactly zero; indeed, one may have Y^' ~ A 4 , Y}i < 0(A 8 ), Y^ < 0(A 4 ), and so 
on (e.g. 0). The solution we give below satisfies these bounds. In writing such forms, it 
is understood that the coefficients multiplying a given power of A may differ from unity, 
but not by as much as a positive or negative integer power of A. The pattern (|29"D is known 



to be experimentally viable [22], and our U(1)a charge assignments constitute a new way 
of obtaining this pattern. 
Consider two cases 

(2a) 9 = -1. Then, as for case (la), 3^ = X y§ + %■ Since Y 2 \ ~ A 4 , Y 2 \ = Y 3 U 2 ~ A 2 , 
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F 3 3 ~ 1 and Y 2 % = Y$ 2 ~ A 4 + < >, we must fit fl29|). We have Y 2 % ~ A 6 + < ^ > and 
we need ~ A 4 , so we require ~ xV (2 ~ b) for & - 0, 1, 2. Since 4a 2 + 2ai = -2, CC 

& + 0'(2-6)-3 = O. (30) 

This cannot be satisfied for 6 = 2, while for b = 1, 9' = 2 which would mean 3^22 ~ x' 1 X 2 or 
Y" 22 ~ A 2 + ... which is too big. Hence, the only acceptable solution is 

6 = 0, &= 3/2. (31) 

We also require Y$ ~ A 6 , i.e. ^i 2 ~ xV (3-<0 for a = 0,1,2,3, and Yf 2 ~ A 6 , i.e. j}f 2 ~ 
X°X^ 3 f° r c = 0, 1, 2, 3. This leads to the CC equations 

a + (3/2) (3 - a) - 3a 2 - 3 = 0, (32) 

c + (3/2) (3 - c) - 3a 2 - 4 = 0, (33) 

which immediately implies a = c + 2. Since a cannot be bigger than 3 we need only consider 
the cases c = or c = 1. Thus we have the solutions 

6 = 0, c = 0, a = 2, a 2 = 1/6, a x = -4/3, 9' = 3/2, (34) 

and 

6 = 0, c = 1, a = 3, a 2 = 0, on = -1, 9' = 3/2. (35) 

The last case allows y™ 2 ~ x' 2 i i- e - ^12 ~ ^ 4 which is too big. Therefore the only viable 
solution corresponds to the first set, eq. (j34]). For this solution, we find 



( A 8 A 6 A 4 \ / A 6 A 4 A 4 

A 6 A 4 A 2 , Y d ~ A 2 A 4 A 2 A 2 I ; (36) 

^ A 4 A 2 1 J \ A 4 A 2 1 



where, for example, the coefficients of Y33, Y 23 , Y" 32 and Y 22 terms could be small enough to 
satisfy experimental bounds, and so forth for other entries. 

Our explanation for property (iii) also requires that the entire contribution to the lepton 
matrices arises from higher- dimension operators. A charge assignment which satisfies this 
requirement is 

Q L + <?e c + <K = 2(ai + a 2 ) - 1, (37) 

a\ = ai = —4/3, a 2 = ct2 = 1/6. (38) 

These equations imply that the charge assignments for the lepton sector are as for the down- 
quark sector, i.e. Y e ~ Y d at Mx- This allows m T pa rrib and mrm^me pa mdm s mb at 
Mx and, given the above-mentioned freedom in the coefficients of the powers of A, this can 
produce a viable model for lepton masses. | 

Given the charge assignments (|38|), the three linear equations (|24]),(|25|) and (p7[) can be 
used to further restrict the solutions to the anomaly constraints (p3|). For example, for the 
first 3-parameter family of solutions we require x = —(l/6)(3z + 7) and y = (1/6) (92 — 13). 

2 Since a factor of 3 falls within the accepted range A < a%j < A -1 , the model of Ref. j23| is an illustration 
of how one can fit experiment with Y e ~ Y d . 
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To pursue this line of research further, the next step is to investigate how the U(1) J 4 
charge assignments that we have made can be derived from a deeper theory (presumably the 
underlying string theory). Another topic for study, but one with much weaker constraints, is 
that of neutrino masses and mixing. Further details will be given in Ref. [20|. This research 
was partially supported by the NSF Grant PHY-93-09888. 
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